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I. INTRODUCTION 


The device-application potential of ultrathin periodically-layered 
semiconductors, commonly known at? super 1 at t ices , has extracted a great 
deal of recent attention, since the pioneering work o i' Esaki and 
collaborators 1,2 in the early 1070's. In constrast to the mainstream of 
research on III — V materials, our interest in semiconductor super iatt icei 
has focused on the possibilities offered by creating such structures in 
silicon, still the most technologically important electronic material. 

In this report we summarize the second phase of a comprehensive 
theoretical study of si 1 icon-superlattice electronic structure and its 
implications for enhanced carrier mobilities in such systems. The 
first, preliminary phase of this work was performed under NASA grant 
NAG-1-33 and included detailed studies of hypothetical one- and 
two-dimensional systems, as well as the formulation of the basic 
three-dimensional model upon which the present investigation has 
f ocused . 

Our efforts during the present grant period have been concentrated 
in five main areas : 

(i) A highly-refined and intensive investigation of both the 
valence- and conduction-band-edge electronic structure of the thin-layer 
(< 11 A) s i 1 icon-super 1 at t ice systems, utilizing the methodology 
developed during the first phase of this work. The results of this 
study have been recently published. 3 

(ii) The development of a new theoretical method for extending our 
thin-layer calculations to layers of thicknesses > 11 A, where most 
potential experimental interest lies. 

(iii) A corresponding intensive investigation of the electronic 
structure of thicker-layer (11-110 A) silicon super 1 at t i ces . 

( iv) Preliminary calculations of impur i ty-scat ter ing-1 imi ted 
electron mobility in the thicker-layer super 1 at t ices . 

(v) The selection and testing of an experimental method for 
producing the fine metal lines that would be required to produce an MOS 
super lattice. 

In addition to the principal investigator, S. Kr i shnamur thy , as a 
graduate student in the Department of Physics, participated in all five 
of the above research areas. The work reported in areas (ii), (iii), 
and ( iv) has comprized a major portion of his Ph.D. thesis 4 and will 
also be published. 6 Additional par t ic ipatants in area (v) were 
Dr. K. P. Roenker, Assistant Professor of Electrical Engineering, and 
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E, Craig, graduate student In the Departmant of Electrical and Computer 
Engineering . 

In Sec, II, the major theoretical considerations which have guided 
our work are reviewed. Then in Secs. Ill and IV, respectively, our 
electronic structure calculations covering areas (i), (ii), and ( l i l ) 
and our mobility calculations covering area (iv) are presented and 
discussed. Finally, our efforts to experimentally produce the fine 
metal lines of area (v) are summarized in Sec. V. Concluding remarks 
are given in Sec, VI, while a bibliography of conference papers and 
journal publications resulting from this work are collected in Sec, VII. 
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II. GENERAL THEORETICAL CONSIDERATIONS 
A. Possible physical structures 

* 

Within the confines of silicon materials technology (i.e., 
excluding 1 1 1-V and 1 1 —V I materials), the number of possible 
superlattice structures is greatly limited. With regard to conventional 
layered structures, the Si-Si ,„ x Ge x system is the prototype silicon 
super 1 at t ice . In addition to our theoretical work on this system, a 
preliminary experimental investigation of Si-Si,_ x Ge x super 1 at t ices has 
been made at Rockwell International. 8 Quite interestingly, their 
results showed evidence of enhanced transverse electron mobility even 
without deliberate modulation doping. This observation is actually in 
sharp constrast to the GaAs-Ga t _ X A 1 x As system, in which the uniformly 
doped superlattice shows significantly decreased mobility below the 
bulk, 7 and strongly suggests that additional physical mechanisms are at 
work in Si-Si j_ x Ge x . Our theoretical studies indeed suggest a second 
possible mechanism, one of quantum origin arising from the uniquely 
modified electronic structure in silicon-based super 1 at t ices . In 
particular, we find a reduced transverse conductivity effective mass 
associated with the superlattice band structure, as discussed in 
Sec. Ill below. It is unclear at this point, however, whether or not 
this effect is actually playing a role in the observed enhanced 
mobility. For one thing, the experimental program at Rockwell 
concentrated on generally thick-layer (> 300 A) systems, where such a 
quantum effect should be greatly diminished. Secondly, the Rockwell 
Si-Si j_ x Ge x super 1 at t i ces have been grown with a chemical -vapor- 
deposition (CVD) technique in which the quality of the interfaces is not 
rigidly controlled. Possible interface strain fields and other 
nonuniformities cloud any straight-forward interpretation of their 
results in terms of our theoretical predictions. Attempts to grow 
better super 1 att ices , utilizing the more-precise technique of molecular- 
beam-exp i taxy (MBE) 3 and accomodating the lattice mismatch in a 
strained-layer structure, 8 have recently been undertaken at Bell 
Laboratories . 9 

Beyond the Si-Sij_ x Ge x system, silicon materials technology offers 
some other novel possibilities for superlattice structures. One is an 
MOS (metal-oxide-silicon) configuration with a grided metal electrode 
made up of fine parallel metal lines uniformly spaced on the oxide 
layer. 10 (See Fig. 1.) In this system, it is envisioned that an extra 
periodic electron potential could be extended into the silicon by 
holding alternate metal lines in the grid at different voltages. The 
primary technological challenge of creating such a system comes in 
producing a periodic metal electrode pattern of the ultrathin dimensions 
required (< 500 A line spacings and thicknesses). We have done some 
exploratory work in this direction, focusing on ( i ) existing techniques 
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Fig. i. Conceptual drawing of a possible MOS superlattice structure, 




for producing fine metal lines of the required thickness and (il) 
methods which could be tested with our laboratory facilities. After a 
literature search of current fine-line technology, we selected one 
promising method for laboratory testing. The results of this study are 
reported in Sec. V. 

Another possible silicon super lattice structure could be one 
consisting of alternately n- and p-doped silicon layers. This type of 
system (the so-called nipi super 1 at t i ce 1 1 ) has been realized in GaAs 
with MBE, 12 and shows a number of interesting properties including a 
tunable band gap. 13 It turns out, however, that advantages such as 
modulation doping in enhancing carrier mobilities are lost in this type 
of super 1 at t ice . 


B, Enhanced mobility 

As in bulk semiconductors, there are three basic physical 
mechanisms which influence carrier mobilities in the ideal semiconductor 
super latt ice : (i) impurity scattering, (ii) phonon scattering, and 

(iii) band edge electronic structure. The first, impurity scattering* 
is primarily dictated by the level and uniformity of doping. Its 
degrading effect on carrier mobility in the superlattice can be 
minimized by modulation doping. In this technique, only one of the two 
constituent layers in the superlattice is doped, and this is done in 
such a way as to segregate the added carriers from their parent ions. 

For example, if one dopes the constituent possessing the highest-lying 
bulk conduction-band edge (e.g., Ga 1 _ J[ Al x As) with shallow donors, the 
ionized donor electrons will migrate to the second constituent (e.g., 
GaAs) having the lower conduction-band edge. Application of a 
transverse electric field (i.e., one parallel to the layers) then allows 
the donor electrons to move in regions free of ionized impurity 
scattering centers, leading to enhanced transverse mobility. This 
effect is particularly striking at low temperatures, 7 where impurity 
scattering dominates t other scattering mechanisms. Furthermore, the 
effect is not dependent on any quantum-mechanical coherence of the 
super lat t ice , but only on the heterostructure itself, and. has now been 
observed in a single-period GaAs-Gaj^Al x As heterostructure as well, 14 

For the Si-Sij^Gejj super latt ice , modulation doping would appear to 
be 'ess advantageous in principle because in this case it is the pure 
material constituent (Si) which posseses the higher-lying 
conduction-band edge, so that now the segregated donor electrons must 
move in the (disordered) alloy component (Sij_ x Ge ). Thus reduced 
impurity scattering is traded for increased disorder scattering. This 
balance may nonetheless still be favorable, however, because the latter 
is not expected to be large, but this has not as yet been demonstrated 
exper imentai I y . 
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The precise role of phonon scattering in altering superlattice 
mobility is also not well established at this point Ln time. Model 
calculation of the electron-phonon scattering rate in a layered 
superlattice simulating GaAs-Gaj.jjAl x As suggests, in fact, that this 
rate should actually be enhanced rather than diminished. 15 The 
calculation assumes among other things, however, that both the phonons 
and the electronic band structure are unaltered in the super latt ice , 
which is, generally speaking, not the case. 

Carrier mobilities are also directly affected by the electronic 
structure of the material in which the carriers move. Of course, part 
of this dependence is ah’eady implicit in the impurity and phonon 
scattering rates, but we know from simple treatments of bulk 
semiconductors that a major effect of the band structure ultimately 
comes in the form of a multiplicative factor involving the conductivity 
effective mass m* , Ideally, for ionized Impurity scattering the 
mobility varies as (m*)" 1 / 2 and for acoustic phonon scattering as 
(m*)" 6 / 2 . In the limit of a three-dimensional quantum-mechanical ly 
coherent super latt ice , where carrier mean-free paths are long compared 
to the superlattice period s and the band structure is well defined, one 
expects similar premiums on effective mass, with m^ for the bulk 
replaced by an appropriate value for the super l at t ice . In silicon 
super 1 at t i ces , this should be the case when s << 1000 A at room 
temperature or s << 4000 A at liquid nitrogen temperature. 

Interestingly, one can make a number of general statements about 
the expected behavior of effective masses in the superlattice prior to 
doing detailed calculations. First, the usual isotropic conductivity in 
bulk diamond or zincblende crystals becomes anisotropic in the 
superlattice with distinct values parallel (i,e., transverse) and 
perpendicular (i.e., longitudinal) to the layers of the super 1 at t i ce . A 
scalar conductivity effective mass in the bulk thus becomes a 
conductivity effective-mass tensor in the super lat t ice , involving both 
transverse and longitudinal masses, say, m^ T and m|J L for holes and rn* T 
and m* L for electrons. As in the bulk, band edges in the superlattice 
occur along symmetry directions, and the electronic structure in the 
vicinity of each edge can be described in terms of transverse and 
longitudinal band masses, say, m^ j and mj| t for the valence band and m* [ 
and m* t for the conduction band. (IJere, as in customary bulk 
descriptions, longitudinal l refers to the direction of the symmetry 
axis and transverse t to the plane perpendicular to -that axis.) In the 
usual case, where there is a simple one-to-one correspondence of the 
band edge in the super lattice with the band edge in the bulk, it is also 
useful to relate these latter masses to their counterparts in the bulk. 
For example, for the heavy-hole valence band with a band edge at the 
center of the Brillouin zone (Jc — 0, the T point), one can write 
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m h l = «| m hh t 1 ) 

and 

m ht = a t m hh » ^ 2 ) 


where m^ h is the bulk heavy-hole mass and a t and a t are enhancement or 
reduction factors which give the relative hole masses in the 
super 1 att ice . Furthermore, in this case 


m hL 


m ht 


(3) 


and 

m hT = m ht • 


Similarly, for a single conduction band edge {not necessarily at ]< = 0), 

m e I = 0 1 m | (5) 


m*t = m* ' < 6 ) 


where mj and m* are bulk mosses and /S ( and /S t are relative electron-mass 
factors. If the conduction band edge is also at It * 0 (i.e., a 
direct-gap material), then mj = m* and also m* L = m* t and m* T - m* t . 

The more important case for our purposes is that of a silicon-like band 
structure with multiple equivalent band edges centered out along the 
<100> axes near the X-point zone boundaries. Then m* T and m* L are 
weighted averages of m* t and m * { . In the bulk, or course, there are six 
equivalent, effectively degenerate, minima with m* = m* T = m* L and 


“• [(l/mp+t 8/m; >] ' 


(7) 
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In the super latt ice , symmetry alone will determine whether or not this 
six-fold degeneracy is lifted and a different average of m* t and m* ( is 
operative. For super latt ices oriented in a <lll> direction, ail six 
minima remain equivalent by symmetry, so that ‘no lifting of the 
degeneracy is expected. For super latt ices oriented in either a <110> or 
a <100> direction, on the other hand, a four-two splitting of the 
six-fold degeneracy is expected. For example, in a [I00]-or iented 
silicon superlattice the [100] and [TOO] minima are equivalent by 
symmetry as are the remaining four minima located in the (100) plane. 
Thus the superlattice is expected to have either two or four equivalent 
minima. Our calculations, discussed in Sec. Ill below, show, in fact, 
that either case can occur, with the two-minima case preferred as the 
thickness of the super latt ice increases. In the two-minima case, one 
has simply 

m e L = m el (8) 


and 

mJ T = m’ t (9) 


in place of an average like Eq. (7). 

One can further anticipate certain general trends in the factors 
a t , atj , /8 t , and /3 j . When the super latt ice is formed, the most 
pronounced effects on the band, structure are the folding of the bulk 
bands into the mlni-Br i 1 louin zone of the super lattice and the opening 
of band gaps (the so-called minigaps) on the longitudinal zone-boundary 
faces. The latter has the clear effect of flattening the bands in the 
longitudinal direction and hence one expects a, > 1 and, for conduction 
band edges 'located along the longitudinal direction, > l. In 
practice, this mechanism is so effective that except for the very 
thin-layer cases one usually approaches the conditions 

Qtj -* °° and jSj -» oo . (10) 


Consequently, the longitudinal conductivity effective masses m^ L and m* L 
are always increased over the bulk, and no mobility enhancement is 
likely. In any transverse direction, on the other hand, such flattening 
of the bands does not occur and one expects to a first approximation a t 
* 1 and /S t 1, so than definite conclusions about m^ t and m* t can't be 
drawn. Interestingly, however, one can easily see that even if ( 3 1 = 1, 
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m* T based on Eq. (9) is significantly decreased below m* from Eq, (7). 
Using the known buik silicon values m] = 0.97m and = 0.19m, 16 one 
f inds m* - 0.20m and 


7t 



= 0.73 . 


( 11 ) 


Thus symmetry alone could account for a 27% decrease in the transverse 
electron conductivity effective mass in [l00]-or isnted silicon 
super 1 at t i ces . Ultimately, this effect arises because of a more 
favorable averaging of the relatively small transverse band mass m| in 
the superlattice than in the bulk. Moreover, it is an effect clearly 
associated with the indirect nature of the bulk band gap and is absent 
in super lattices derived from direct-gap materials such as GaAs . 

Of course, additional considerations enter the picture beyond 
band-edge effective masses. For example, the bulk conduction-band edges 
which split away and above in energy in the super lattice will always be 
accessible at a high enough temperature. In practice, a thermal average 
over a number of bands is necessary to establish the effective carrier 
mobility. Moreover, a different thermal average, which weights the low 
energy states more heavily in the superlattice than in the bulk, is 
usually appropriate. In the bulk, the conduction-band-edge 
constant-energy surfaces are ellipsoidal in shape and the density of 
states for a single band has the familiar form 


N C (E) 



m;(2m;) l / 2 (E-E c ) 1 / 2 


( 12 ) 


where E c is the conduction-band edge. In the super 1 at t i ce , on the other 
hand, the band flattening mentioned above makes the constant-energy 
surfaces more nearly cylindrical in shape, so that the corresponding 
density of states is proportional to that of a two-dimensional electron 
gas. For a single conduction-band edge at J< = 0: 


h 2 


rrs 


N C (E) = 


0t m * • 


(13) 


which is independent of energy. Our preliminary calculations of 
impur i ty-scatter lng-I imi tod electron mobility in silicon superlattices 
based on these considerations will be discussed In Sec. IV below, 

♦ 

C. Band-structure model 

In formulating a band-strucutre model for silicon super lattices, we 
have been guided by a desire to examine general trends in this new class 
of systems rather than to obtain highly accurate details in any specific 
structures. Towards this end, we have adopted an empirical 
tight-binding (TB) approach which both gives a reliable bulk band 
structure of pure silicon and provides enough flexibility to treat all 
superlattice structures of interest with the variation of only a few 
parameters. The specific super lattice geometry considered here Is 
illustrated in Fig. 2. We assume that the underlying tetrahedral 
coordination of pure bulk silicon is preserved in the super latt ice , so 
that the position and bond length d of the nearest neighbors is 
unchanged. The shaded areas in Fig. 2 represent regions where either 
the silicon atoms have been subst i tut ional ly replaced by a second 
species (i.e., germanium in Si 1-x Ge x ) or the electron potential has been 
rigidly raised or lowered as envisaged in MOS or nip i superlattice 
structures. In the former case, we further ignor any disorder in the 
alloy by making a virtual crystal approximation and also any slight 
Increase in the bond length due to the presence of the second species. 
Then in all cases the periodicity of the superlattice is unchanged in 
the transverse plane (the yz plane in Fig. 2), but in the longi t iudinai 
direction the layers repeat with a period 


s = n fl a/2 , 


(14) 


where a is the bulk silicon lattice constant and 2n a is the number of 
atomic layers in each primitive cell of the super 1 att ice . For 
convenience, we assume n s to be an even integer > 2, Then both the 
superlattice and its minizone in reciprocal space have simple tetragonal 
symmetry, as shown in Fig. 2. In the TB method, a small number of 
atomic-like basis functions is assigned to each atom and the electronic 
states of the solid are obtained as linear combinations of these 
orbitals. Following Harrison, 17 and the earlier work of Chadi and 
Cohen 18 on bulk silicon, we have adopted a minimal basis set of one 
s-like and three p-1 ike orthonormal atomic orbitals. A careful study of 
the electronic structure of pure silicon revealed that in this basis one 
must retain both nearest-neighbor and second-neighbor interatomic 
interactions to obtain a satisfactory description of the energy bands in 
the vicinity of the fundamental energy-band gap. The near-neighbor 
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Superlattice Geometry 



Fig. 2. Schematic drawing of the silicon superlattice geometry treated 
in this work. At the top is the real-space structure with 
period s; below ic the corresponding m i n i -Br i 1 1 ou i n zone 
(minizone) in reciprocal space. 
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interatomic matrix elements In covalent solids are primarily {unct .ums 
of bond length d In the solid and, according to Harrison, 17 may btf 
expressed in the form 


V U'm = ^U'm h 2 /(md 2 ) 


{16) 


where 1,1' * $ or p and the subscript m is either a or rr. We have fixed 
the dimensionless parameters Vii’n together with the Intra-atomic 
energies e g and e p by fitting our bulk TB energy bands to the accurate 
pseudopotential results of Chellkowsky and Cohen 10 for energies up to 1 
eV above the conduction-band edge. The best set of values we have 
thereby obtained are listed in Table l of Ref. 3. 

Since we assume the bond length d of the superlattice to be 
identical to that of bulk silicon, the V^, m are also unchanged In the 
super 1 att ice , The silicon intra-atomic energies £ a and e p , on the other 
hef>v become 

£ s + V< a > ( 16) 


and 

tp + V<P> (17) 


for the second layer of the super 1 att ice , Fitting of the bulk silicon 
energy bands fixes £ a and e p , so that the entire effect of the 
superlattice potential is incorporated into the constants and V^P). 

In general, these quantities represent different potentials seen by s 
and p electrons in the second layer. The special case 

v s = V< 8 > = V(P> ( 18) 


corresponds to the situation envisaged in an MOS or nipi super latt ice , 
where the electron potential in alternating layers is raised or lowered 
by a constant value V s , We have designated this as type-I behavior. 3 
Another interesting case occurs for 


V = - v< 3 ) = P * 
o s 8 * 


( 19 ) 
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in which case s and p electrons effectively see equal and opposite 
superlattice potentials. Detailed consideration of the electronic 
structure of bulk germanium and the measured band-edge discontinuities 
of the Sl-Ge heterostructure 20,21 suggests that this limit, which we 
have designated as type-II behavior, 3 is appropriate to the Si-Si l-x Ge x 
super I at t ice . Our studies have further shown that in the absence of 
band bending at the heterojunction interface that an approximate 
relationship between the superlattice potential V s and the germanium 
concentration is 


V a * (0.5 ± 0. l)x (in eV). 


( 20 ) 


Band bending and other nonldeai interface effects may be incorporated 
into our model by allowing, for example, V s to grow linearly across 
several layers at the interface instead of changing abruptly from zero. 
We have considered a number of such graded-layer calculations, but they 
produce results so qualitatively and quantitatively similar to the ideal 
case that we shall not consider this complication further in the present 
di scuss 1 on. 
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III. ELECTRONIC-STRUCTURE CALCULATIONS 

In the above energy-band model, the superlattice electronic 
structure is reduced to a function of two chemical parameters, Vj' 3 ) and 
V<P) , and three geometrical quantities; superlattice orientation, 
period s, and relative thickness of the two alternating layers. We have 
focused our attention on the limiting cases of highest applied interest: 
a [ I00]-or i ented superlattice of either the type-I or type-I l variety 
with equal layer thicknesses. The latter condition is maintained to 
maximize the effective superiattice strength for fixed values of V s and 
n s . Unequal layer thicknesses produce qualitatively similar but 
quantitatively smaller effects in all cases. Also note that for equal 
layer thicknesses l, 

s = 2 l , (21) 


and the parameter n g is just the number of {100}- atomic planes contained 
in a given layer. 


A. Thin-layer super! att ices {l < 11 A) 

In a straight-forward application of the TB method, the energy 
bands E(jt) are obtained by direct diagonal izat ion of the appropriate 
Hamiltonaln matrix. In either the bulk or the superiattice, the 
Hamiltonian matrix is constructed from a basis of Bloch wavef unct ions of 
the type 


*!<<*<£) = vFf? 


( 22 ) 


where a is an orthonormal s, p x , p y , or p z atomic orbital and the sum is 
over the N equivalent sites R t on which this orbital is centered in the 
crystal. If there are 2n a atoms in the primativp cell, then there are 
8n„ such basis states and an Bn„ x 8n„ Hamiltonian matrix results. The 
detailed form of this matrix is given and discussed in the Appendix of 
Ref. 3 . 


We have performed TB electronic structure calculations for the 
superiattice band energies E(J<) with k in both a longitudinal [100] and 
a transverse [OOtJ direction in type-I as well as type-I I cases over the 
parameter ranges 
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o < v < V? ax 

8 8 


( 23 ) 


and 


2 < n g < 8 , 


(24) 


where V£ aX is the value of V s for which the fundamental band gap closes. 
For n a > 8, matrices larger than 64 x 64 would result and alternate 
means of calculation using complex band structures are needed, as 
discussed in Sec. 1 1 1 B below. The limit n s = 8 corresponds to layer 
thicknesses of Z = 11 A or a superlattice period of s = 22 A. 

Typical results that we have obtained for the band structure in the 
vicinity of the fundamental energy band gap are illustrated in Fig. 3. 
The top panel shows the pure bulk silicon bands with the characteristic 
indirect gap and conduction-band edges along the <100> directions out 
toward the zone boundaries. The super lattice band structure may be 
thought of as evolving from these bulk bands in two steps, The first 
step is the J$-space translation or folding back of the bulk band 
structure into the minizone of the superlattice for infinitesimal V s . 

The second is the additional dispersion introduced into the folded bands 
for finite V g . The center panel of Fig. 3 shows the bulk silicon bands 
folded into the n s = 4 superlattice minizone. Four of the 
conduction-band minima remain fixed, but the remaining two along the 
[100] and [TOO] directions are displaced to positions near & = 0. When 
the superlattice potential V g is increased from zero, the bands split 
apart, the fundamental gap narrows, and the [100] and [TOO] minima are 
further displaced in position and also lowered in energy relative to the 
four transverse minima, as shown in the bottom panel of Fig. 3. 

The valence-band edge in the superlattice remains stationary at 
= 0, the minizone center, except for very large values of V s (e.g., V g > 
2.0 eV for n g = 2). The conduction-band edge J< c , on the other hand, 
varies with both n„ and V„, with either the characteristic 
two-longitudinal-minima or four-transverse-minima behavior occuring in 
each case. The qualitative position of ]c 0 as a function of n a and V g 
for both type-I and type-I I super 1 at t ices is illustrated in Fig. 4. 

Note that the two-minima case is favored ts n s increases and that in 
certain cases the two minima coalesce at - 0 and a direct band gap 
results. The general trend suggests that i sufficiently large n g , the 

band gap will become direct for all values or V g . Our thick-layer 
calculations confirm this, as discussed below. 
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The magnitude of the fundamental band gap E g decreases continuously 
as V s is increased in all cases. The gap decreases more rapidly with V B 
the larger n s is and also somewhat faster in the type-II case than in 
the type-I case, as shown in Fig. 5 for n s = 2 and n Q = 8. The latter 
behavior is expected on the basis of Eqs, (16 )— ( 19), since the s and p 
energies are pushed together in the second layer for the type-II 
super lattice. 

At the valence- and conduction-band edges ve have also determined 
the longitudinal and transverse curvatures of the energy bands to obtain 
the effective masses m^ j , m^ t , m* j , and m* t . Relative band masses otj , 
a t , /S| , and /? t were then obtained by normalizing the calculated band 
masses by the calculated values of m£ h , mj , and m* for bulk silicon via 

Eqs. (1),(2),(5), and (6). The relative masses so obtained show a 

number of interesting and prominent features, as illustrated in 
Figs, 6-6. The valence-band masses CX| and a t are smooth functions of V s 
for fixed n g since the band edge ic y is stationary, As expected, we find 
a t > 1 for V a > 0, with otj >> 1 possible for large n s , The transverse 
mass a t , on the other hand, drops off very sharply with V s from unity to 
about 0.8 in all cases, as shown in Figs, 6 and 7. For larger values of 
V s , a t remains less than one and is found snaller in the type-II case 
than in the type-I case for fixed n s and is larger the larger n s is for 

both cases. Unlike a ( and a t , the conduction-band masses /? ( and /S t are 

not smooth functions of Vj. becasue j£ c also varies with V g . The expected 
behavior /Sj > 1 and /S t 1 is, however, still found. Figure 8 displays 
/? t versus V a for n s = 2 and n s = 8 in the type-II super lat t ice . 

Whenever the conduction-band edges lie along the longitudinal x 
axis and Eq. (11) applies for the relative transvese electron 
conductivity effective mass 7 T , we indeed find y T < 1, so that enhanced 
mobility is possible. This situation occurs in our calculations for n g 
> 6 in the type-I super lattice and n s > 4 in the type-II super 1 at t ice . 
Figure 9 shows the variation of y T , 7 T “ 1 / 2 , and 7 T " 6/|/3 with V for n s = 

8 in the type-II super 1 at t ice . The quantities 7 T -1 ^ 2 and 7 T ~ 5//2 are the 
ideal enhancement factors expected for impurity-scattering-limited 
electron mobility and phonon-scattering-limited electron mobility, 
respectively, if all other things are equal. 



-18- 


R6H 


I 

i 

r J 

Y_ 


2.0 


a>. 


- 1.0 

2.0 

1.0 

« 

0.0 

- 1.0 


[ 001 ] 


Silicon 


[ 100 ] 


ms 

m 

r t 


£jg 

8 

[001] 

[100] 

fm 

p 

mm 

|^S 

[001] 

[100] 


Superlattice 

n s »4 

Vs-0.0 


Superlattice 
ns“4 
V s *1.0eV 
Type II 


Fig. 3. Evolution of the super lat t ice electronic structure in the 

longitudinal [100] and transverse [ 00 1 ] directions from the 
bulk energy bands for a typical case. Top panel, the bulk 
silicon bands; center panel, folded bulk bands in the n g = 4 
superlattice minizone; bottom panel, type-II superlattice bands 
for a finite superlattice potential V g . 


Wt #* « 






19- 



n a *4 



i 

ii 


n s -e 



V S (©V) 


Fig. 4. Qualitative position of the conduction-band edge in both type-I 
and type-I I silicon super la tt ices as a function of n s and V . 
The position Jc c * [OOlJ actually corresponds to four equivalent 
edges symmetrically centered along the the transverse y and z 
axes; J< c = [100] corresponds to two equivalent edges 
symmetrically centered along the *x axis near \c * 0; = 0 

corresponds to the special case where the latter two edges 
coalesce into a single point at the center of the minizone. 
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Fig. 5. Variation of the fundamental energy band gap E g with V a for n s 
= 2 and n s = 8 in type-I and type-I I super 1 at t ices . The 
vertical arrows indicate places where the conduction-band edge 
position J< c changes di scont inuous ly (see Fig. 4). 
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Fig. 6. Relative transverse hole band mass a t vs V s for n s = 0 in 
type-I and type-I I super 1 at t ices . 
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Relative longitudinal a. and transverse a. hole band masses as 


a function of n and V in the type-II super 1 at t ice . The scale 


at the bottom of the figure relates the concentration x to V 


via Eq. (20) . 







Transverse: n s -8 


x in Si-Si i_ x Ge 


Relative transverse electron conductivity effective mass y T and 
impur i ty— scat ter ing—1 imi ted and phonon— scat ter 1 ng- 1 imi ted 
enhancement factors y T -1 / 2 and 7 T “ 5 / 2 vs V 3 for n s =* 8 in the 
type— 1 1 super latt ice . The lower scale is as in Fig. 7. 
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B. Thicker-layer superlattices (11 £ l 5 110 A) 

Because the TB Hamiltonian matrix Increases In size proportional to 
layer thickness, the calculation of the superiattice electronic 
structure in this manner becomes intractable for even modestly thick 
layers. The complex band-structure technique of Schulman and Chang 22 
nicely overcomes this difficulty. In this method, one views the 
superiattice as being composed of bulk regions In each layer separated 
by interface planes, as shown in Fig. 10, The number of atomic 
Interface planes is taken such that the atoms in bulk regions I and II 
are not directly coupled: for second-neighbor interactions, two planes 

at each Interface are required. One then capitalizes on the fact that 
the solutions to Schrbd inger ' s equation are local and hence in the bulk 
regions must be some linear combination of bulk states inclucing both 
extended Bloch states (real Jfcc) and evanescent states (complex J<) . For a 
given energy B, there are in our case 16 possible J< states In each bulk 
region. These states, together with two-dimensional Bloch sums of 
atomic orbitals centered on the interface-plane atoms, can be used to 
expand the superiattice wavefunction in the form 


+ . (25) 


where q * q x ^ + jcg is the superiattice wavevector, (j * nsjc, with n an 
integer, = k x £ + JC| ( is a bulk wavevector (1 = 1,16), is the 

corresponding bulk wavefunction for region a (cr = I, II in ^lg, 10), and 
a is a two-dimensional Bloch state for interface plane j (j = 
S lt S 2 ,S 3 ,S 4 in Fig. 10) and atomic orbital a (a « s ,p x ,p y ,p z ) . The 
resulting Hamiltonian matrix is of size 48 x 48 independent of layer 
thickness . 

The price that must be paid for a fixed Hami 1 tonian-matr ix size is 
an energy dependent basis, requiring an iteration procedure to calculate 
the superiattice band structure se 1 f-cons i stent iy . For given values of 
n s and V s , one fixes the superiattice wavevector q and selects a trial 
energy E = E t to to be used in Eq . (25). For the chosen E t , the 16 bulk 
states for each region are obtained. This information is contained in 
the complex band structure of the bulk material, 23,24 which need only be 
generated once for each layer. Figure 11 shows the [100] complex band 
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Fig, 10. Schematic diagram of the n a * 0 superlattice geometry, showing 
the projection of four planes of atoms onto the z * 0 plane. 
The unhatched circles represent the Si atoms of layer A, while 
the hatched circles represent the atoms of layer B. Atoms 
belonging to the planes z = -a/4, a/4, and a/2 are marked 
+ , and ++, respectively. The bulk regions I and II and 
interface layers S t , S 2 , S 3 , and S 4 apply to the complex 
band-structure technique discussed in the text. 
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structure of bulk silicon. The 48 x 48 superlattice Hamiltonian Is thon 
set up and diagonalized to obtain the true superlattice energy E B for 
the assumed basis. One Iterates this procedure until E a « E t to some 
acceptable toleronce. In our procedure, an accuracy of 10“ 5 eV In 
energy is achieved with typically 5 Iterations in less than 1 minute of 
cpu time, By varying q, the complete superlattice hand structure may be 
mapped out in this way. 

Our tight-binding complex band-structure method (TB-CBM) Improves 
upon that of Schulman and Chang 22 In two important ways. First, our 
tight-binding model includes second-neighbor interactions as well as the 
first-neighbor interactions considered by them. Second, in constructing 
the superlattice Hamiltonian from the basis (25), we have included 
certain finite-layer corrections 4 ' 5 to the bulk regions that are 
necessary t,o make the method exact and were neglected by Schulman and 
Chang, These corrections become vanishingly small as n s -» <», but 
significantly affect fine details of the band structrure such as 
effective masses for the range of n a considered here. 

TB-CBM calculations of the silicon superlattice band structure have 
been carried out for 0 < n s < 80 ( 11 < l < 1 10 A) and 0 < V s < t ,0 eV 
for both the type-I and type-1 1 cases. In particular, the magnitude of 
the fundamental band gap, conduction band-edge positions, and band-edge 
effective masses have been obtained over these parameter ranges. For n a 
= 8, the results obtained are Identical to those from the thin-layer TB 
calculations, confirming the exactness of our approach. For n a > 8, the 
trends extablished in the thii.-'^yer calculations are continued. The 
variation of the band gap magnitude E with V a and n a is shown in 
Figs. 12 and 13 for thicx-layer lype-T and type— 1 1 super iatt ices , 
respectively. As expected, E g continues to decrease with Increasing V a 
and n a , but the E g versus V a curve quickly saturates for large n s . 
Extrapolation of the n B = 80 curves in Figs. 12 and 13 suggests that the 
band gap will close for V a = 1.11 eV, the bulk silicon band gap, in the 
limit n a -♦ <», as required in the type-I case. 

The positions of the two n s — 8 conductions band edges remain along 
the x axis for n a > 8. For n s £ 32 ( l > 44 A), the two edges 
permanently coalesce at P and a stable direct band gap results. This 
situation is illustrated in Fig. 14, where the complete longitudinal 
[100] and transverse [001] band structure E(jq) in the vicinity of the 
fundamental gap is plotted in a typical thick-layer case. The energy 
difference A between the X and P conduction band edges always remains 
small, however, (5 0.0065 eV = 75 K) and eventually declines back to 
zero as n g °o, as shown in Fig. 15. 

For n a > 8, the infinite limits of Eq. (10) have been essentially 
reached for the longitudinal electron and hole effective masses at any 
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f ini t e V g . The transverse effective mass factors a t and ft t also 
saturate i) values near unity for n 8 > 64. The thick-layer, type-I and 
type-II superlattice values of a t vs V a are plotted in Fig. 16, while 
the thick-layer, type-1 1 values of ft t vs V s are plotted in Fig. 17. The 
cort esponding type-1 values of ft t are virtually Identical to the latter. 
With ft t * 1, one has y T 'v 0.73 and the impurity-scattering-limited and 
phonon-scattering-limited enhancement factors, and y T ~ 5 ^ 2 , 

saturate at values of 1,17 and 2.20, respectively, as shown in Fig, 18. 
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Fig. 11. Complex band structure of bulk silicon in the [100] direction 
with |c = k x £. Left panel: complex bands with Re(k x ) = 0; 

center panel: real bands, with Im(k x ) = 0: right panel: 

complex bands, with Re(k x ) = 2ir/a. 
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Fig. 14. Longitudinal [100] and transverse [001] band structures near 
the fundamental band gap of a type-I I super lattice with n s = 
64 and V. = 0.2 f eV. The bands shown are those thermally 
accessible to mobile charge carriers at room temperature. 
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V, <«V) 

Fig. 15. Energy difference A between the X and T-point conduction-band 
minima vs V s for various type-I and type-II super lat t ices . 
(Solid lines, type-I and dashed lines, type-II.) The n a = 64, 
type-II result has been omitted for clarity. 
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Impur i ty-scatter ing-1 imi ted and phonon-scattering-limited 
enhancement factors 7 T ”*/ 2 and 7 T “ 5 ^ 2 vs V s for various n s in 
the thick-layer, type-I I super 1 at t ice . The n s = 16 result for 
y T ~ 1 / 2 has been omitted for clarity. Compare with Pig. 9. 
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IV. IMPURITY-SCATTERING-LIMITED ELECTRON MOBILITY 

In our present studies, we have considered ionized impurity 
scattering of electrons from the screened Coulomb potential 


v | ( r ) = - 



(26) 


where k is the dieietrlc constant of the host and X" 1 is the Debye 
length, X 2 = (4rrn,e 2 /A:k B T) , for uniform doping to a concentration n t . 

In bulk silicon, where the density of states (12) Is operative for each 
conduction band edge, the corresponding mobility of nondegenerate 
electrons is well approximated by the standard form 26 


32 k z (k B T) 3 / 2 
" V8rr a e 3 


Gg ( n j , T ,m c ) , 


(27) 


where G B is a slowly varying function of nj , T, and m* . In the 
super 1 at t ice , a similar calculation, but with a constant density of 
states like (13), yields for the i-^A conduction band 


k Q T 


. 15S IC 3 / Z «Vni „ . _ , . 

Msl = 4rrV2e 2 (n^*®) 1 / 2 SL<ni ' ’ 


(28) 


neglecting interband scattering. Here G SL is also a slowly varying 
function and m* T is the appropriate transverse conductivity mass for the 
jib band. The inverse-square-root effective-mass dependence of the bulk 
is recovered, as expected, but the temperature and impurity 
concentration dependencies now reflect the two-dimensional nature of the 
transport. Summing contributions from all thermally accessible 
conduction bands, the average electron drift mobility in the 
superlattice is given by 





n i^S L./ n I ’ 


(29) 
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4 

where n, la the electron concentration In the band. For T < 300 K, 
the conduction bands which must be included in Eq. (39) are shown in 
Fig. 14. 


For a typical silicon superlattico, the calculated relative 
mobility based on Eqs . (27)-(29) is plotted in Fig. 19 as a 

function of nj for various values of T and in Fig. 20 as a function of T 
for various values of n ( . It can be seen that at low temperatures, 
where impurity scattering dominates other scattering mechanisms, an 
enhancement in electron mobility is indeed possible for sufficiently 
high carrier concentrations. The qualitative variations found with n t 
and T are most easily understood by considering the high temperature 
limits of Eqs. (27) and (28), from which one finds 



where Cj and C 2 are constants. For a fixed T in the range of interest, 
this function increases to a maximum for some n, and then decreases 
monotoni cal ly . Similarly, for a fixed nj, this function increases to a 
maximum for some T and then decreases. 

The corrresponding Hall mobility in the superlattice may be 
calculated by replacing Eq. (29) by the average 

2 ^^(/^l) 8 /2 n,/4 L , (31) 


where r t is a weighting factor reflecting the energy dependence and 
anisotropy of the scattering. Our calculated relative Hall mobility is 
plotted in Fig. 21 as a function of nj for various values of T and in 
Fig. 22 as a function of T for various values of n If for the same 
superlattice parameters as in Fig. 19. Similar magnitudes and trends 
are clearly seen. 

The r°sults presented in Figs. 19-21 should be considered as upper 
limits to what could be expected in real systems. Neglected factors 
will tend to reduce the actual mobility, in particular, interband 
scattering at high carrier concentrations and alloy scattering in the 
case of a Si-Si j_ x Ge x super 1 at t ice . At the same time, these results 
take only partial advantage of the unique reduced-effective-mass 


39 - 


enhancement mechanism in silicon super 1 att ices , which should be much 
more effective in the case of phonon scattering. 




Fig. 19. Variation of the impur i ty-scat ter ing-1 imi ted electron drift 
mobility with donor concentration n t at various temperatures 
for a type-II super lat t ice , where n a = 64 and V s * 0.2 eV, 
relative to bulk silicon. 
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Fig. 80. Variation of the impur i ty-scat ter ing-1 imi ted electron drift 

mobility with temperature T for various donor concentrations, 
as in Fig. 19. 
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Fig. 21. Variation of the impurity-scattering-limited electron Hall 

mobility with donor concentration nj at various temperatures 
for a type-II super lat t ice , where n a = 64 and V s = 0.2 eV, 
relative to bulk silicon. 
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as in .«* lg . 21. various donor concentrations, 
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V. THIN METAL LINE FABRICATION 

We hu7e made a reasonably thorough literature search of recent 
developments in fine-line technology and have uncovered a wide range of 
techniques all of which can basically be catagorized as either X-ray, 
electron-beam, ion-beam, or optical lithographies. In our laboratory, 
only optical lithography was available, and the most promising optical 
method we discovered was one reported by Speidell 27 at IBM. He claimed 
lines as narrow as 800 A with his technique, and, as a preliminary test, 
we decided to try it in our lab. We closely followed his procedure, 
modifying steps only according to the facilities available to us. The 
complete procedure and results are described as follows. 

On several clean silicon wafers, oxides of thickness *3 3400 A were 
grown, Standard clean 1 {SCI) was prepared as a 7:3:3 solution by 
volume of deionized (Dl) water, H g 0 2 (30%), and NH 4 0H (29%). This 
solution was qsed at about 80°C for 20 minutes to clean the oxidized 
wafers. Next the wafers were prebaked at 180°C for 30 minutes, coated 
with hexame thy 1 di s i 1 azane at 4000 rpm for 15 seconds, then coated with 
AZ 1350J photoresist at 5000 rpm for 20 seconds, and finally baked at 
65°C for 25 minutes. Since we desired an undercut in the photoresist 
layers, the soft-baked wafers were next soaked for 5 minutes in 
chlorobenzene and soft-baked again for 30 minutes at 65°C. 

With a previously-prepared mask for 2.5 /xm lines, the wafers were 
UV exposed for 15 seconds. We used high-speed developer (1:3.5 of AZ 
developer and water) to develop the wafers. They were developed for 9 
sec with 3 minutes of DI rinse, then placed in the developer for 1 
minute and rinsed ‘in running DI water for 3 minutes. The dried wafers 
were hard-baked for 30 minutes at 125°C in vacuum. 

Figure 23 shows SEM photos of typical photoresist structures we so 
achieved. Because of the chlorobenzene soak, undercuts or "lips" are 
clearly visible on either side of the line. These "lips" are less than 
1/3 of’ the total thickness of the line. We estimate that each undercut 
is approximately 4000 A wide. The size of the undercut can be reduced 
by reducing the soaking time in chlorobenzene. The undercut 'is made at 
a rate of approximately 800 A per minute of soaking time in 
chlorobenzene. 

The desired thin metal lines are achieved by evaporating Ai at an 
angle of 20° into the undercuts. To etch away the unwanted Al and 
photoresist, we tried two different methods. First, the wafers were 
heated in acetone for 5 minutes and cleaned in an ultrasonic cleaner for 
l minute. The photoresist surface is thereby removed without disturbing 
the Al . Then an etch was done to remove Al layers between the two lines 
that are formed in the two undercuts. The etch removed Al between the 
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lines but also smeared out the two lines at the same time. To get 
sharper lines, we followed the reverse procedure, That Is, we first did 
the A1 etching and then removed the photoresist. The latter gave 
reasonably sharp A1 lines. 

The final line thickness we achieved was ** 4000 A, dictated by the 
size of the undercuts. This dimension can definitely be reduced to the 
the order of 600 A, as demonstrated by Speidell 27 , by appropriately 
reducing the soak time in chlororbenzene , To obtain even thinner lines, 
on the order of 100-200 A, however, will undoubtedly require X-ray or 
electron-beam lithographies. For this reason and because we' had no 
immediate scheme in mind to produce the periodic line pattern ( p lg. 1) 
required for an M0S super latt ice , we did not pursue this study further. 



VI. CONCLUSIONS 


Our results bode well for the propspect of achieving enhanced 
electron mobilities in silicon superlattice systems, We have seen that 
the general trends in band structure indicate a reduced electron 
conductivity effective mass in the transverse direction, with a unique 
symmetry lowering in [l00]-or iented super lattices playing a leading 
role. For a fixed superlattice potential strength, this effect 
initially increases with layer thickness, but quickly saturates to a 
limiting value, Phonon-scattering-limited mobility, which is dominant 
at room temperature, should be most strongly enhanced by this mechanism, 
with an ideal limiting enhancement of 2.20 predicted if all other things 
are equal. The two-dimensional nature of the transport, however, also 
brings with it a strong dependence of mobility on carrier concentration. 
Our calculations show that impurity-scattering-limited mobility, which 
is dominant at low temperatures, will actually be diminished at low 
carrier concentrations and enhanced only for high carrier 
concentrat ions . 

Since our band-structure model for silicon superlattices is an 
idealization we expect that any mobility calculated on the basis of our 
model should represent an approximate upper limit for given geometrical 
conditions. Two principal effects which'peed to be addressed more 
thoroughly in any refined theoretical treatment are (i) chemical 
disorder in the Sij_ x Ge x layer of the superlattice in the type-I I case 
and (ii) surface effects, including band-bending, strain fields, and 
possible defects at the interfaces between layers. Physically, both of 
these effects should decrease any measured electron mobility. In 
addition, interband scattering, which was neglected in our 
impurity-scattering-limited mobility calculation should also decrease 
the actual mobility. 

The most pressing need for the near term, however, is some hard 
experimental data on high-quality superlattice structures. Hopefully, 
this will be forthcoming in the Si-Si 1 _ x Ge x case with the new work 
initiated at Bell Laboratories . 9 Also, the fabrication of an MOS 
super lattice should be given some serious consideration in light of the 
continuing rapid developments in fine-line technology. 
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